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Is dark energy an effect of averaging? 
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The present standard model of cosmology states that the known particles carry only a tiny fraction of total 
mass and energy of the Universe. Rather, unknown dark matter and dark energy are the dominant contribu- 
tions to the cosmic energy budget. We review the logic that leads to the postulated dark energy and present 
an alternative point of view, in which the puzzle may be solved by properly taking into account the influence 
of cosmic structures on global observables. We illustrate the effect of averaging on the measurement of the 
Hubble constant. 



1 Introduction 

The concordance model of cosmology states that only 4% of the mass/energy content of the Universe 
is carried by identified particles, and the dominant contributions are unknown dark matter (22%) and 
mysterious dark energy (74%) ||T|. It is said that we are Uving in the age of "precision cosmology", which 
is justified as a statement on the precision of observations, but certainly not on their understandings. The 
purpose of this talk is to recapitulate the genesis of the dark energy hypothesis and to show that it poses a 
serious crisis to theoretical cosmology. We argue that a possible resolution for this crisis may come from 
the proper treatment of effects from averaging. 

2 Dark energy 

The vast majority of cosmic background photons (photons without identified source) are due to the cosmic 
microwave background (CMB) radiation, which is isotropic (apart from a component due to our proper 
motion) with fluctuations at the level of 10~^. Thus the isotropy of the Universe (around us) is an obser- 
vational fact. On the other hand, there is no evidence that the Milky Way would be a distinguished place 
in the Universe and thus it seems reasonable to follow Copernicus to assume that we are not living at the 
centre of the Universe. Therefore, we end up with a homogeneous and isotropic model. The most general 
line element in such a Friedmann-Lemaitre model is given by 



where a = a{t) denotes the scale factor and K/a^ the spatial curvature {K = —1, 0, +1) of the Universe. 
In this model physical distances evolve with time, rph = af{r) (we may fix the present scale factor 
ao = 1), where /(r) = {Arsinh(r), r, arcsin(r)} for K = {—1,0, 1}. It is useful to define the rate of 
linear expansion H = a /a. Without using Einstein's equations, just assuming that photons travel on null 
geodesies, one finds that light from distant objects is red-shifted by z = AA/A = 1/a — 1 and the so-called 
luminosity distance (defined by the observed flux from a source with known luminosity c?l = \/L/ 4ttF) 
satisfies the Hubble law 
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Fig. 1 The binned magnitudes A/i (see text) for different SN data sets and light-curve fitting methods. 



at small redshifts z ^ 1. Observations show that Hq ~ 70 km/s/Mpc > and thus we know that the 
Universe expands. The typical time scale of the cosmos is estimated as in ~ 14 Gyr, in concordance 
with the estimate on the age of the oldest objects in the Universe [21. Another interesting kinematic 
quantity is the acceleration (deceleration) of the cosmic expansion, encoded by the deceleration parameter 
q = —(a/ a) / H'^ . Measurements of supernovae of type la (SN la) indicate that < 0. 

One might wonder what could be the reason for the high symmetry of the observed Universe. The 
answer of modern cosmology is the idea of cosmological inflation [31, an epoch in the very early Universe 
that is held responsible for the Universe's homogeneity and isotropy. On top of that, generic predictions of 
cosmological inflation are spatial flatness (JIk = 0) and the existence of tiny fluctuations ID, as observed 
at a level of 10"^ in the CMB. 

We may ask how strong is the evidence for the acceleration of the Universe, even before we wish to 
invoke a model for the matter content, without assuming the validity of Einstein's equations. This test goes 
as follows. We ask at what level of confidence we can reject the null hypothesis that the Universe never 
accelerated (at z < few). For a spatially flat model, as suggested by the idea of cosmological inflation, our 
null hypothesis can be expressed as 

dL < ^(l + z)ln(l + z). 
-no 

Based on two different SN la data sets, two different fitting methods and two different calibration methods, 
we conclude that the null hypothesis is rejected at > Str 131. 

This test can also be done in a calibration-free way, i.e. without assuming certain values for the Hubble 
constant Hq and the luminosity of the SNe. Instead, we consider the averaged value of the apparent 
magnitude m of SNe within a certain redshift bin relative to the averaged apparent magnitude mnoarby of 
nearby SNe with z < 0.2. Defining a new quantity 

A/.t = m - m„oarby = 5 logio "^L - 5 logio dL.noarby, 

the null hypothesis can be written as A/i < 0. For the test, we average m within redshift bins of width 
0.2. In Fig. [U the result is shown for two data sets (Gold Q and ESSENCE l7j) and two fitting methods 
(MLCS2k2 m and SALT ||9l). The evidence for acceleration can be obtained by dividing A/i by its error. 
The result for each redshift bin is given in Table[T] Also given is the evidence obtained from averaging over 
the apparent magnitude of all SNe with redshift z > 0.2. It is noticeable that using the calibration-free test 
weakens the evidence for acceleration (4.3(t) compared to the result of the previous test (5.2(7). 
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Table 1 Statistical evidence for acceleration from different redshift bins for different SN data sets and light-curve 
fitting methods. 

Let us now ask what that implies for the dynamical model of the Universe. The most general equation of 
motion for a homogeneous and isotropic Universe under the assumption of diffeomorphism invariance and 
at most second order derivatives on the metric (Lovelock theorem lilOJ ) provides us with the Friedmann 
equation and the continuity equation (c = 1) 

7?2 + ^-A=^6, e + 3i/(e+p) = 0. 
3 

e and p denote the energy density and pressure. G and A are Newton's constant and the cosmological 
constant. In order to close this system of equations, we also need to impose an equation of state p = p(e), 
which in many cases can be expressed as w = p/e. A dimensionless energy density 51 = e/(3i?^/87rG') is 
introduced for convenience. 

We have already mentioned that cosmological inflation predicts a vanishing spatial curvature for all 
practical purposes and thus we can put K = 0. In the present Universe, matter is non-relativistic (cold), 
which means p <C e and thus we put p = 0. If we now would put A = as well, we arrive at the Einstein- 
de Sitter Universe, which predicts q(t) = 1/2 and to ~ 2tH/3 ~ 9 Gyr. Therefore, we cannot stick to this 
simple model, as we have seen that q < at high confidence level and the oldest stars are known to be as 
old as 12 Gyr or even older IS). 

These facts became clear only about ten years ago — the 1998 cosmology revolution. One can rewrite 
the dynamic equations into 

-3- = 47rG(e + 3p) - A < 0. 
a 

The inequality and thus acceleration, can occur if either A > 47rGe (equivalent to JIa > rim/2) or there 
exists some form of energy density with p < — e/3, a form of energy that has never been observed in the 
laboratory so far and is thus dubbed the "dark energy". For p = — e this reproduces the behaviour of the 
cosmological constant (as e = in that case), which is thus a special case of dark energy. 

So a minimal model consistent with the predictions of cosmological inflation and allowing for the 
observed accelerated expansion of the Universe is given by A' = p = and JIa > flm/2 = (1 — riA)/2 
and thus JIa > 1/3 and n,n < 2/3. This model is called the flat ACDM (cold dark matter) model or 
concordance model, as it is currently in good agreement with all cosmological observations. 

One can test for a possible deviation from A' = or p = by means of CMB observations, SN la and 
the baryon acoustic oscillations and finds that there is no observational evidence to question the validity of 
the flat ACDM model I.Lillj. 

Let us stress that an important assumption in the reasoning described above is the homogeneity and 
isotropy of the Universe and thus it is also important to test this assumption with the help of SN la. A 
simple test can be done by fitting the Hubble diagram (c^l vs. z) for SNe from just one hemisphere and 
comparing the result of the fit to the opposite half of the sky lfT2l . Weinhorst and one of the authors 
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have analysed four different SN data sets and find a statistically significant deviation from the isotropy of 
the Hubble diagram. However, it seems that these anisotropies are aligned with the equatorial coordinate 
system (aligned with the Earth's axis of rotation) and are thus Ukely to reflect systematic problems in 
the SN la observation or analysis. Another recent analysis revealed anisotropies in the Hubble diagram 
obtained by the HST Key Project [1131 . If these anisotropies would be confirmed by future SN surveys and 
a systematic effect could be excluded, this would pose a serious challenge to the interpretation of Hubble 
diagrams in the Friedmann-Lemaitre context. 



3 Problems of the concordance model 

Thus it seems that the concordance model is able to describe everything, however there are three conceptual 
problems in this model. 

The cosmological constant problem. We miss out a fundamental understanding of the vacuum energy 
of quantum fields. This vacuum energy density is indistinguishable from a cosmological constant and thus 
we have introduced an element to the Universe that we do not understand. The solution for this issue must 
be expected from a future theory of quantum gravity, as gravity (i.e. geometry of space-time) is the only 
"force" that couples to the energy density of the (quantum field) vacuum. Thus this problem is obviously 
beyond the scope of what we can hope to understand currently. 

The coincidence probleni(s). The second problem however is a pure cosmological one. It is called 
the coincidence problem and there are two ways to formulate it. The first formulation has some anthropic 
touch, but is widely used: how does it come that 17,„(to) ^ i^A(io)? This is indeed surprising. Assume we 
would have lived when the Universe had half (twice) of its present linear size [at z = 1 (-1/2)], the ratio 
^m/^A would be larger (smaller) by almost an order of magnitude. It is therefore quite surprising that we 
observe f^A ~ 0.7 > 1/3 at to- 
One can reformulate the above situation via replacing mankind by large-scale-structures on the typical 
scale for structure formation. According to inflationary cosmology, the primordial fluctuations in matter 
and metric are free of any scale. The only scale that enters in the physics of the large scale structure 
is the horizon scale at matter-radiation equality. This scale appears because the radiation pressure acts 
against the gravitational instability in the early Universe and structure growth starts only once the pressure 
support fades away. Thus the shape of the power spectrum of density fluctuations is changing at this scale. 
According to the concordance model, the horizon scale at matter-radiation equality is seen at the 100 Mpc 
scale today. Scales above that show the primordial slope of the spectrum, scales below that exhibit a 
substantially decreased slope. It is now interesting to observe that those structures grow to typical density 
contrasts of order {Se/e)z=o.x ~ (1 + Zcq){6e/e)eq ~ 0.1 and it turns out that z = 0(0.1) is the epoch 
when the typical scale of hierarchical structure formation starts to evolve non-linear (i.e. it decouples from 
the overall cosmological expansion). The coincidence is that this is also the epoch when acceleration seems 
to set in 1141 . Note that these large structures are indeed observed, e.g. the Sloan great wall and various 
big voids, as big as a few 100 Mpc ITSl . 

The averaging problem. This is an old question in general relativity lfT6l . Solving the Einstein equa- 
tions for an averaged space-time is certainly not the same as averaging the Einstein equations for a very 
inhomogeneous and anisotropic space-time, because these two operations, i.e. averaging and time evolu- 
tion do not commute. Thus the issue is, how big is this effect and what are its observational signatures. 
Again the scale of interest must be the 100 Mpc scale, as on much smaller scales, say 1 Mpc objects like 
galaxies have already presumably undergone the process of violent dynamical relaxation and formed sta- 
tionary bound systems that cannot be treated like dust particles anyway. However it is clear that for objects 
like the Sloan great wall, which just started to collapse yesterday, there might be important effects showing 
up in the averaging process. 

In summary, these three aspects are entangled with the issues of inhomogeneities and anisotropies in 
the local Universe. The beginning of the domination of dark energy coincides with the onset of structure 
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formation. Consequently, light may be shed on the dark energy crisis (i.e. the accelerated expansion of the 
Universe) by studying the averaging problem in the perturbed space-time. 

An attempt for setting up an formaUsm to handle this problem has been pioneered by Buchert ifTTlfTSl . 



4 Cosmological backreaction 

These problems provide a good motivation to consider the possibility that dark energy is not a fundamental 
component of the Universe, but our averaged description of the Universe is not appropriate for a correct 
interpretation of observations lfT4l [191 l20l . Indeed many cosmological observations are averages. Two 
important examples are the power spectrum P{k), which is a Fourier transform and thus a volume average 
weighted by a factor cxp(ik • x), and the Hubble constant Hq. Let us consider an idealised measurement 
of Hq II2TI . One picks N standard candles in a local physical volume V (e.g. SN la in Milky Way's 
neighbourhood out to ~ 100 Mpc), measures their luminosity distances di and recession velocities Vi = 
czi, and takes the average Hq = J2iLi 'T- limit of a very big sample, it turns into a volume 

average Hq = y J f d^. For objects at 2; ^ 1, the spatial average is appropriate for the average over 
the past light cone, as the expansion rate of the Universe is not changing significantly at time scales much 
shorter than the Hubble time. 

We now recapitulate the formulation of Buchert, which we will use below to argue that the averaging 
effect is indeed sizeable and does give important contributions to cosmology: 

Buchert's setup is well adapted to the situation of a real observer On large scales, a real observer is 
comoving with matter, uses her own clock, and regards space to be time-orthogonal. These conditions 
are the definition of the comoving synchronous coordinate system. Buchert used physically comoving 
boundaries, which is the most natural approach in this setup. In the following the Universe is assumed to 
be irrotational as a consequence from cosmological inflation. 

In synchronous coordinates, the metric of the inhomogeneous and anisotropic Universe is ds^ = —dP + 
gij {t, x.)dx'^dx^ , and the spatial average of an observable 0{t, x) in a physically comoving domain D at 
time t is defined as fTTl 



(0>D = -^ / 0(t,x)v/detg~dx, (1) 



D 



where Voit) = Jjj y^detg^dx is the volume of the comoving domain D. We may introduce an effective 
scale factor qd iflTl 



/ T/ \ 1/3 



aoo \ Vdo 



(2) 



The effective Hubble expansion rate is thus defined as Hd = do/aD = {())d/3 (9 being the volume 
expansion rate). 

From the definition ([T]), we obtain effective Friedmann equations ifTTI from averaging Einstein's equa- 
tions, 

doY SttG Hd inG , o ^ 

foff, = (ecff + 3pcff), (3) 



^a£) ) 3 au 3 

where t^s and pcft are the energy density and pressure of an effective fluid, 

ecft^(e>D-Y^((Q)D + (i?)i5), V.^^-^:^[pD-\{R)i^. (4) 

((5)u = — (6')|,) — 2(0-^) u is the kinematical backreaction (cP' being the shear scalar), and (i?) d 

the averaged spatial curvature. They are related by the integrability condition ifTTll 

{a%{Q)D) +a%(a\{R)u) (5) 
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We define the equation of state for the effective fluid as Wcs = Pcs/^cB- It is highly remarkable that any 
spatially averaged dust model can be described by an effective Friedman-Lemaitre model. 

We can map this effective fluid on a model with dust and "dark energy". Let n be the number density 
of dust particles, and m be their mass. For any comoving domain (n) o = {n) Do {(IDo/o,d)^ ■ For a dust 
Universe, in which e{t, x) = mn{t, x), we identify em = {e) o — m{n) o, and from ^ the dark energy is 
hence described by edo = d + {R)d)/{^Q'^G). From ^ we find that constant (Q) d = — (-fi')_D/3 

corresponds to the case of a cosmological constant A = {Q)d- Equations (O and (01 are not closed and 
additional input is required. Below we close them by means of cosmological perturbation theory. 

To study the scale dependence of physical observables {Q)d, {R) d, (e) d, Hd and Wcb, we calculate 
them to second order in a perturbative series of the effective scale factor qd- We start from a spatially flat 
dust model. Its scale factor a{t) is different from the effective scale factor a^i in (O, and their relation was 
provided in ll22ll . In synchronous gauge, the Unear perturbed metric is 



where 'J and x are the scalar metric perturbations, Dij = didj — ^SijA and A denotes the Laplace 
operator in three-dimensional Euclidean space. The solutions for 5' and x are given in terms of the time 
independent peculiar gravitational potential </j(x): = ^AtptQ^^t'^^^ + |(^ and x = Stptf/'^t'^^^ (only 
growing modes are considered) 1221 |20l . ip is related to the hypersurface-invariant variable C |[23j by 

Following II22I . with the help of the integrability condition, we yield the scale dependence of the aver- 
aged physical observables up to second order 1241 {{0)di = Jjj Odx/ dx hereafter) 

{Q)d = ^B{^)tl, (6) 
an 

{R)d = ^^{A^)Di-5^B{^)tl (7) 
3 ajj an 

{^)d = T^^, (8) 



Hn = 



r, 3/2 



(9) 



WeS = -—4{Aip)m-^4{Biip)--{A^)j,,\, (10) 

with B{ip) = {d'{ditpAip) - d'idjipd^d^ip)) di - |(Av3)|,^. We see from Q - ^ that these quantities 
are functions of surface terms only, so all their information is encoded on the boundaries. The dependence 
on cosmic time of these averaged quantities can be found in 1221, and their leading terms are gauge invariant 



We now turn to the estimate of the effects of cosmological backreaction as a function of the averaging 

1/3 

scale r VJ^^ . We show below that cosmological averaging produces reliable and important modifi- 
cations to local physical observables and determine the averaging scale at which 10% corrections show 
up. 

Effective acceleration of the averaged Universe occurs if Ccff + 3pcff < 0, i.e. (Q) o > 47rG(e) d- Thus 
we estimate 



{Q)i 



47rG(e) 



^ "ii Tj,.^.4 8 Rfi 8 1 fRa^* 



2^^(^)^" = 27 (TTIF^^^^ ^ ^(TT^ I V j '''' 



i?H — 3.00 X lO'^h ^ Mpc is the present Hubble distance and T'^ = 2.35 x 10 ^ the dimensionless power 
spectrum |[T| (C ~ —5lp/3 on superhorizon scales). In ( fTTT i. we use a = 1/(1 + 2) and to = 2i?H/3, 
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since B is already a second order term. Each derivative in B is estimated as a factor of 1/r in front of the 
power spectrum. As Lp is constant in time, we identify today's Vtp with 97^^/25 at superhorizon scales. If 
I (Q) D /47rG(e) d\ > 0.1, cosmological backreaction must be taken into account seriously. From (fTTT i. we 
find 



21h-^ Mpc 



For observations at z <C 1, < 30 Mpc {h — 0.7). 

The criterion for the scale at which the effect of the averaged spatial curvature (i?) £> emerges is deter- 
mined by 



Eeff 



167rG(e)z3 



1 



3 1 



r 



We expect a 10% effect from (i?) d within 



rR < 



At small redshifts, r/? < 77 Mpc. Moreover, a 1% effect is expected up to scales of ^ 240 Mpc. Note 
that the curvature of the Universe has been measured at the few per cent accuracy in the CMB [I]. It was 
shown in ll26l that even small curvatures might influence the analysis of high-z SNe significantly. 



5 Signature of cosmological backreaction: the effective Hubble rate 

It is now interesting to discuss the fluctuation of the local measurement of the Hubble expansion rate. 



Hd — Ho 



2 



The leading contribution to this effect involves the averaged curvature terms and therefore cannot be thor- 
oughly described in a Newtonian setup 1271. An effect larger than 10% shows up for 

VI + z 

which reads rn < 54 Mpc at small redshifts. 

Let us compare (T% with observations from the HST Key Project ||281 . We use 68 individual mea- 
surements of Hq in the CMB rest frame from SN la, the Tully-Fisher and fundamental plane relations 
(Tables 6, 7 and 9 in 1281 ). The respective objects lie at distances between 14 Mpc and 467 Mpc. As 
(fT2] | can be trusted only above 30 Mpc, we drop the 4 nearest objects. Be r^. Hi and (Ti the distance, 
Hubble rate and its la error for the i'th datum, with increasing distance. We calculate the averaged dis- 
tance for the nearest k objects by ffc = 9i^i/ 5*' ^i'^^ weights gi = l/crf . An analogue holds 
for the averaged Hubble rate Hk, i e. Ho for different subsets. The empirical variance of each subset is 
tr^ = J2i=i 9iiHi — Hk)'^ /[{k— 9i]- We stress that (fT2] | is insensitive to global calibration issues. 

The comparison of our result (fTZt with the HST Key Project data is shown in Fig. 2. The "global" 
Hq is the central HST value 72 km/s/Mpc 1281 . We find that the theoretical curve (oc l/r^) matches the 
experimental data fairly well, without any fit parameter in the panel. Since ( fT2] i is just an estimate, we vary 
the result by factors of 0.5 and 2, which gives lower and upper estimates to the theoretical curve. In Fig. 2, 
the experimental data and their error bars are in agreement with these estimates. 

Before we can claim that we observe the expected 1 /r^ behaviour of (fTSI) and thus evidence for cos- 
mological backreaction, we have to make sure that the statistical noise cannot account for it. The expected 
noise is given by means of the variance of ( fT2b in the situation of a perfectly homogeneous and isotropic 
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Fig. 2 The scale dependence of the normalised difference between the effective Hubble rate Hd and its "global" 
value Hq = 72 km/s/Mpc (Data from the HST Key Project (28)). We expect that all data points should lie around the 
two thick lines (~ ±l/r ) given in l ll2t . Thin lines indicate the theoretical error of our estimate. Dashed lines are the 
expected statistical noise for a perfectly homogeneous and isotropic model. 



model with measurement uncertainties <Ji. In the case of a perfectly homogenous coverage of the averaged 
domain with standard candles, we expect a l/r^/"^ behaviour. In Fig. 2, we show the statistical noise for 
the actual data set, which turns out to be well below our estimate ( fTSl i and the data points. At ^ 50 (80) 
Mpc, the value of Hjj differs from its global one by about 10% (5%), whereas the expected variance for a 
homogeneous and isotropic Universe is only 3% (1%). 



6 Toward a new model and open issues 

What can we conclude so far? We showed that the large scale structure affects the observed Hubble 
expansion rate and that the averaged spatial curvature leads to sizeable corrections up to scales of about 
200 Mpc. Consequently, we must expect that also other cosmological probes at z < 1, like galaxy clusters 
or galaxy redshift surveys, could be affected by cosmic averaging. 

Is that enough to explain dark energy? Currently, this is not the case. In our framework an effective 
acceleration is possible below scales of 20 Mpc. But SN la observations show strong evidence for the 
acceleration of the Universe, based on objects much more distant than that. It seems to take more to 
explain dark energy. A quite attractive idea is a Universe in which large voids dominate the Universe's 
volume, see 129 , 30l [3n[32l . These works rely on well motivated toy models. It remains to be shown that 
these void dominated cosmologies follow from the inflationary paradigm. 

On the other hand, our analysis is based on an almost scale-invariant power spectrum predicted by 
inflationary cosmology, but is limited to z ^ 1, as the average must be taken on the past light cone and not 
on a spatial slice (see, e.g. 133)1 ). A formalism (like Buchert's one for spatial averaging) for the light-cone 
averaging of an arbitrary inhomogeneous and anisotropic space-time is needed before we can know if dark 
energy is merely an illusion due to averaging. 
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